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ABSTRACT 



Aims. We use the parallax measurements of RE J 0317-853 to determine its mass, radius, and cooling age and thereby constrain its 

evolutionary origins. 

Methods. We observed RE J 0317-853 with the the Hubble Space Telescope's Fine Guidance System to measure the parallax of 

RE J 0317-853 and its binary companion, the non-magnetic white dwarf LB 9802. In addition, we acquired spectra of comparison 

stars with the Boiler & Chivens spectrograph of the SMARTS telescope to correct the parallax zero point. For the corrected parallax, 

we determine the radius, mass, and the cooling age with the help of evolutionary models from the literature. 

Results. The properties of RE J 0317-853 are constrained using the parallax information. We discuss the different cases of the core 

composition and the uncertain effective temperature. We confirm that RE J 0317-853 is close to the Chandrasekhar's mass limit in all 

cases and almost as old as its companion LB 9802. 

Conclusions. The precise evolutionary history of REJ 0317-853 depends on our knowledge of its effective temperature. It is possible 

that it had a single star progenitor possible if we assume that the effective temperature is at the cooler end of the possible range from 

30 000 to 50 000 K; if T^jf is instead at the hotter end, a binary-merger scenario for RE J 03 17-853 becomes more plausible. 
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1. Introduction 

REJ 0317-853 is a unique hydrogen-rich white dwarf, which 
was discovered as an EUV source by the ROSAT Wide Field 
Camera (|Barstow et al.||1995|l. An analysis of follow-up spec- 
troscopy established that the stellar surface is covered by a very 
strong magnetic field with a range of about 170-660 MG, imply- 
ing that REJ 0317-853 has one of the strongest magnetic fields 
detected so far in a white dwarf 

The optical spectrum together with UV observations taken 
with the lUE satellite and the Hubble Space Telescope indicated 
that REJ 0317-853 possesses a very high effective temperatures 
in the range from 30000 to 55 000 K; [Barstow et al] ( [T995] l 
achieved their best fit for about 49 000 K. A careful analysis of 
the EUVE spectrum using the interstellar medium Lyman lines 
to account for the interstellar extreme ultraviolet absorption im- 
plied an effective temperature of 33 800 K ( [Vennes et al.|2003] l. 
Within these constraints, REJ 0317-853 is one of the hottest 
known magnetic white dwarf (MWD); in any case, it has the 



highest known temperature of all MWDs with a field strength 
above 20 MG ( jKawka et al.|2007l|Kiilebi et al.|2009| l. 



Barstow et al. ( 1995 1 performed high-speed photometry 



demonstrating that the optical brightness of REJ 0317-853 
varies almost sinusoidally with a period of 725.4 ±0.9 sec and 
an amplitude of more than (JVl; these results were confirmed by 
Vennes etal. (2003), who inferred a period of 725. 727 ±0.001 sec 
from the variation in the circular polarisation. The only rea- 
sonable explanation of these results is rotation, implying that 
REJ 0317-853 rotates more rapidly than any other known white 
dwarf that is not a member of a close binary. The photomet- 
ric variation must be caused by differences in the brightness on 
various parts of the stellar surface. Since no strong absorption 
lines are detected in the optical, a possible explanation may be a 
variation in the effective temperature over the stellar surface; the 
reason for this temperature inhomogeneity is currently not well 
understood but is probably connected to stronger or weaker con- 
tributions to the magnetic pressures in the stellar atmosphere at 
different locations on the stellar surface with different magnetic 
field strengths. 



* Based on observations made with the NASA/ESA Hubble Space 
Telescope, obtained at the Space Telescope Science Institute, which is 
operated by the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS5-26555. The Guide Stai' Catalogue- 
II is a joint project of the Space Telescope Science Institute and the 
Osservatorio Astronomico di Torino. 



a clearer insight into the evolution of 
Burleigh et al. ( 1999[l obtained phase- 



To achieve 

REJ 0317-853, 

resolved far-UV Hubble Space Telescope (HST) Faint Object 
Spectrograph spectra. They found that the previous optical 
results could generally be confirmed, but that the splitting of 
the Lyman a component into subcomponents implied that the 
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field is probably more complicated than indicated by the mean 
optical spectrum. By compiling a time series of spectra, a model 
for the magnetic field morphology across the stellar surface 
was produced using the radiation-transfer models through a 
magnetised stellar atmosphere from Jordan (1992] see paper 
for a basic description of the methods) and an automatic least 
squares procedure. The magnetic geometry could be equally 
well described by an off'set magnetic dipole (xoff - 0.057, 
yofs = 0.004, and Zo« - -0.220 stellar radii), which produces a 
surface field strength distribution in the range 140-730 MG or 
an expansion into spherical harmonics up to / = 3 in which the 
surface field strengths are constrained to be within the range 
1 80-800 MG. 

The mass of the white dwarf was constrained by estimating 
the absolute magnitudes (or absolute fluxes) calculated from the 
spectroscopic fit parameters T^ff, log g and white dwarf evolu- 
tionary models (e.g. Wood 1995; Be nvenuto & Althaus||1999] l. 
The determination of the mass of RE J 0317-853 is, in general 
not straightforward because of the effects of strong magnetic 
fields; the usual method of using the Stark broadening of the 
spectral lines to determine log g and subsequently a mass-radius 
relation fails in the presence of a magnetic field of several hun- 
dred MG; the reason is that the standard theory for Stark broad- 
ening assumes degenerate energy levels but the magnetic fields 
help remove this degeneracy. 

Nevertheless, the mass determination procedure of 
REJ 0317-853 can be improved by the knowledge of its 
distance. REJ 0317-853 is inferred to be in a wide-binary 
double-degenerate system due from its visual companion, which 
is a non-magnetic DA white dwarf companion (LB 9802) 7" 



Table 1. Spectroscopically derived parameters of LB 9802. 



away. This object was analysed initially by Barstow et al. 
( |1995|l, t hen later by Kawka et al. ( 2007| i (for fit parameters see 



Tablelljl. Barstow et al. ( 1995 1 derived a distance in the range 
33-37 pc with these parameters using the evolutionary models 
( 1992 1. The physical companionship of LB 9802 and 



of Wood 



REJ 0317-853 has recently been confirmed by Spitzer IRAC 
obsevations ( |Farihi et al.||2008| ) that demonstrated the common 
proper motion nature of the system. 

With an effective temperature of 50 000 K and assuming a 
distance of 36 pc Barstow et al. ( 1995 i concluded that the ra- 
dius of RE J 0317-853 is about 0.0035 R0 with a corresponding 
extrem e mass of 1.35 M0(logg - 9.5). Later Vennes & Kawka 
([2008^ derived a mass of 1 .32 + 0.03 M© using T^ff = 33 800 K, 
iogg = 9.4 and 27 pc for the distance. If these conclusions are 
true, REJ 0317-853 would not only be one of the hottest known 
MWD but also the most massive (^ 1.35 Mq) isolated (due 
to the large separation of REJ 0317-853 and LB 9802 we can 
assume that both stars did not interact during stellar evolution) 
white dwarf discovered so far; only two other white dwarfs are 
known with masses in excess of 1.3 Mq: LHS 4033 with a mass 
in the range 1.31-1.34 Mpl jDahn et al.|2004"| l an d the magnetic 
white dwarf PG 1658+441 with L3 1 + 0.02 M© (jSchmidt et al" 



[1992 ). 

From the theory of stellar evolution, there are two different 
ways to produce these massive white dwarfs: either by single- 
star evolution of a star with an initial mass higher than 7 or 8 M© 
( [Dobb ie et al. 2006; Casewell et al. 2009, Salar is et al |2009l l or 
from the merger of two white dwarfs with C/O cores (see e.g. 



[Segretain et al. 



rapid rotation o 



1997 1. The latter scenario is supported by the 



REJ 0317-853. 



Jordan & Burleigh' (1999) measured the circular polarisation 



to have a degree of 20% at a wavelength of 5760 A, the strongest 
ever found in a MWD. Together with the assumed small ra- 
dius and strong gravity in the stellar photosphere, this also made 



Ref. 


V 
/mag 


/K 


logs 


dL 
/pc 


1 

2 

3 


14.11 
13.90 


16 030 ±230 
16 360 ±80 
15 580 ±200 


8. 19 ±0.05 
8.41 ±0.02 
8.36 ±0.05 


33-37 

30 

27 



lBarstoweta"L] ( |l995t ; ^Ferrario et a"r] ( |l997| l; ^ [Kawkaet aLl ( |2007t 



REJ 0317-853 a test object for setting limits on gravitational 
birefringence predicted by theories of gravitation, which violate 
the Einstein equivalence principle ( Preuss et al.|2005 i. 

Since the mass determination of REJ 0317-853 was based 
entirely on the uncertain spectroscopic distance of the system, 
we applied for observing time with the HST to measure the 
trigonometric parallaxes of the white dwarf binary system to 
either confirm or disregard the conclusions of Barstow et al. 
(1995). In this paper, we present the analysis of the parallax mea- 
surement with HSTs Fine Guidance Sensor (FGS). 



2. Observation 

2.1. Observations with fhe FGS of fhe HST 

The observations of the magnetic white dwarf REJ 0317-853 
(ffiCRS = 03M7™16U750, 5„„/c«s = -85°32'25'.'45) and 
its non-magnetic white dwarf companion LB 9802 (aicRS = 
03''17'"19^.3050, S,.„icRs = -85°32'31'.'15) with the Hubble 
Space Telescope were performed with the Fine Guidance Sensor 
Ir (FGS Ir) at thi-ee epochs (March 2007, September 2007, and 
March 2008, see Table[2]). The Fine Guidance Sensor is a two- 
axis, white-light shearing interferometer that measures the angle 
between a star and HSTs optical axis by transferring the star's 
coUimated and compressed light through a polarising beam split- 



ter and a pair of orthogonal Koesters prisms (see Nelan et al. 
|1998[ |Nelan|[2010| for a description of the instrument design). 
When FGS Ir is operated as a science instrument, HST pointing 
is held fixed and stabihzed by FGS2 and FGS3 which operate as 
guiders. 

To derive an astrometric solution for position, proper mo- 
tion, and parallaxes, REJ 0317-853, LB 9802, and the reference 
field stars had to be observed at a minimum of three epochs, 
preferably at the seasons of maximum parallax factor to allow 
us to cleanly separate their parallaxes from their proper motions. 
These seasons are separated by about six months. 

Fortunately, the epochs of maximum parallax factor also re- 
sulted in HST roll angles (which are constrained by date) such 
that the two white dwarf stars and the optimal set of astrometric 
reference stars could be observed at all epochs. Figurefl] shows 
the parallactic ellipse and the orientations of the FGS aperture 
at the times of the observations (there were two March epochs, 
2007 and 2008). Experience shows that a minimum of two orbits 
per epoch are required to achieve the highest possible accuracy 
in the final parallaxes. Tablel2] provides the dates of the six or- 
bits for our HST programme. Since the two white dwarfs are 
only a: 7" apart, we were able to use the same reference stars 
for the two white dwarfs using no more HST orbits than would 
be necessary for a single parallax measurement. Using identical 
reference stars also ensured that the parallax difference between 
the two putative companion stars was measured more precisely 
than their absolute parallaxes since the measurements share the 
same correction of relative to absolute parallax. In addition, their 
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Table 2. HST orbits for i/5r proposal 10930 and 11300. 



proposal ID 



start time 



end time 



10930 Mar 24 2007 17:54:01 

10930 Mar 24 2007 19:29:48 

10930 Sep 27 2007 03:28:55 

10930 Sep 29 2007 01:47:19 

1 1300 Mar 29 2008 02:00:30 

11300 Mar 29 2008 03:36: 19 



Mar 24 2007 18:53:25 
Mar 24 2007 20:29: 12 
Sep 27 2007 04:28:19 
Sep 29 2007 02:46:43 
Mar 29 2008 02:59:53 
Mar 29 2008 04:35:42 
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( |1998| l. Since the Pickles library does not cover all spectral sub- 

types, interpolation by eye was performed where appropriate. 

visit l"br late G- and especially K stars, the MK class III templates 
were also looked at because in some cases the star actually 
turned out to be a giant; giants can be clearly distinguished from 
dwarfs, which show an indentation at 5200 A that the giants do 
not or only slightly exhibit. The few metal-poor and metal-rich 
templates were also used, although the difference in the Pickles 
spectra is too small to really make a discrimination in this re- 
spect. 

The absolute magnitude determination was based on an in- 
terpolation of the data taken from |Lang| ( | 1 992| l and Allen's astro- 
physical quantities ( Cox 2000 1. It was achieved by parametrising 
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02 




Palpha 



Fig. 1. The parallactic ellipse of the RE J 0317-853 field and the 
orientation of the FGS Ir field of view at the dates of the obser- 
vations. The X-axis of the FGS Ir is nearly parallel to the line 
connecting the circles that mark the epochs at which the obser- 
vations were made. 



relative proper motions can be measured to provide an additional 
check on whether or not the two white dwarfs constitute a bound 
pair. 

2.2. Spectroscopy of the astrometric reference stars 

Since only relative parallaxes can be measured with HST, we had 
to estimate the parallaxes of a sample of reference stars in the 
vicinity of our target objects which comprise our local reference 
frame (see Figure|2]l. Ref4 and Ref5 were not observed by the 
FGS Ir since they were not needed. 

Spectra of these surrounding stars of similar (or some- 
what larger) brightness than RE J 0317-853 were taken in ser- 
vice mode with the Boiler & Chivens spectrograph of the 
1.5m SMARTS telescope, located on Cerro Tololo at the 
Interamerican Observatory in Chile, in two nights between 
February 16 and 18, 2008. To ensure that the whole optical range 
is covered, we performed exposures with two gratings (9/Ic and 
32/Ib). Both observing nights were affected by passing clouds 
and the relatively high airmass (> 1 .8) due to the large declina- 
tion difference between the observatory's zenith and the target 
field. Since this could not fully be corrected by flux standards, 
the energy distribution in the blue channel may be compromised. 

The classifications for the reference stars were performed by 
comparing the flux-calibrated spectra to the templates of |PickIes| 



the spectral type so that spectral type F corresponds to 0, G to 1 , 
K to 2, and M to 3, and the spectral type subdivisions correspond 
to the first decimal, i.e. an G2 star would be represented by 1.2. 
A 5th degree polynomial is then fitted to determine the My - 
spectral type relation shown in Figure|3] This was performed for 
both luminosity class III and V, assuming that all our stars come 
from these two luminosity classes. The absolute magnitudes of 
the reference stars were then calculated using these two func- 
tions with their spectral class parametrised in the same way as 
an argument. 

The determination of the errors is not straightforward, since 
not all error sources can be easily quantified. The error in the de- 
termination of the spectral type was roughly quantified by calcu- 
lating the absolute magnitude of the spectral subtypes closest to 
the determined ones were calculated using the same fit function 
(for those stars where the derived spectral type was in-between 
two subdivisions, i.e. in the cases of reference stars 1, 3, and 
9 the second next subtype was chosen). The difference between 
this absolute magnitude and the absolute magnitude obtained for 
the star is then our estimate of the error in the absolute magni- 
tude caused by the uncertainty in the spectral classification. This 
assumes that the error in the spectral type is not larger than one 
subdivision, which might not be true in all cases but should gen- 
erally be the case. It was generally found that the difference in 
absolute magnitude between the measured spectral type and its 
neighbours is about 0.2 mag, so this value was assumed in all 
subsequent calculations. This error of 0.2 mag corresponds to 
an error of 12% in distance (see Table[4j 7th column). The cor- 
responding error in the parallax was used for the correction of 
the relative parallaxes. Given the relation between parallax and 
distance, the error in the former is not symmetric if that of the 
former is. The asymmetric nature of the parallax error is repre- 
sented in column 8 of TableH] The errors given in TableH] do 
not represent the overall error. A main source of error will most 
likely be the photometry, which is not of the highest precision. 
Moreover, our spectra do not allow us to determine the exact 
evolutionary status of the objects, which influences the accuracy 
of the absolute magnitude. For the same reason, the influence 
of metallicity cannot be taken into account, and all stars are as- 
sumed to be of solar abundance. Adding these uncertainties with 
some margin leads to an overall error in distance of 20-30%, with 
the stars Ref7-9 having the larger errors, since we only have one 
spectrum (red of Ref7, blue for the other two) for these objects. 
Since the parallax is the reciprocal of the distance, the stars with 
a large distance are the more reliable ones, especially the two 
giants (Ref3 and 8). 

3. Analysis of the FGS data 

Our astrometric measurements used FGS Ir in position mode 
to observe RE J 0317-853, LB 9802, and the associated refer- 
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Table 3. Coordinates and photometry of RE J 0317-853, LB 9802, and the reference stars. 



name 




Q'ICRS 


^ICRS 


V 


B- V 


U-B V-R V-1 


GSC2 F^ 










/mag 


/mag 


/mag /mag /mag 


mag 


RE J 0317-853 


03" 17"" 16"; 1750 


-85°32'25'.'45 


14.90 ± 0.02 


-0.16' 


-1.13' 0.01' -0.11' 


15.09 


LB 9802 




03''17™19!3050 


-85°32'31'.'15 


14.11 ±0.02 


+0.07' 


-0.68' -0.06' -0.18' 


14.22 


RefP 




03''20"1Z918 


-85°34'56'.'175 


9.42 


0.38 






Ref2'' 




03''18'''52'-;01 


-85°35'20'.'8 


12.27 


0.50* 




12.55 


Ref3 




03''18"'03U 


-85°36'02" 


14.60 


1.14^ 






Ref6 




03''13"'59';7 


-85°30'16" 


14.00 


1.04^ 






Ref7 




osMs-^ssg 


-85°30'20" 


15.00 


0.84' 






Ref8 




0346""4ff.3 


-85°29'48" 


14.37 


1.10^ 






Ref9 




OS^lS-'SSl 


-85°36'42" 


14.36 


0.63^ 
















=HD 2329 
MK class, 






' From'Barstow et al. 


11995 


, ^'http://tdc-www 
etical B -V values 


liarvard.edu/catalogs/gsc2.html ; ^ 
interpolated for spectral type and 


8=TYC9495-788-li|H0g et al. 


1998 ; 


^ =GSC0949500756; 


^ theor 


see Tablelil 





Table 4. Results of the analysis of the reference stars {V magnitude, spectral type and luminosity class, absolute magnitude, distance 
modulus, distance, quantifiable distance error (see text), and parallax and its errors (the subscripts "-" and "H-" refer to the distances). 



star 


V 


spectral 


Mv 


my - Mv 


B-V 


d 


6d 


'^tl 




/mag 


type 


/mag 


/mag 


/mag 


/pc 


/pc 


/mas 


Refl 
Ref2 


9.95 

12.27 


F3-4V 
F7V 


3.48 
3.95 


6.47 
8.32 


0.41 
0.50 


197 
461 


24 
55 


5-08^055 
2j7+o:» 

1 5+^-^^ 


ReO 


14.60 


K1-2III 


0.48 


14.12 


1.14 


6668 


800 


Ref6 


14.00 


K4V 


6.96 


7.04 


1.04 


256 


31 


■i Q 1+0.53 

]1+^-^^ 


Ref7 
Ref8 
Ref9 


15.00 
14.37 
14.36 


KOV 

Kim 

G3-4V 


5.98 
0.55 
4.81 


9.02 

13.82 

9.55 


0.84 
1.10 
0.63 


637 

5808 

813 


76 

720 

98 



Table 5. Astrometric results for RE J 0317-853, LB 9802. and the astrometric reference stars: parallax, proper motion in right 
ascension, proper motion in declination, the standard errors of the proper motion, and the standard errors in the fiducial coordinates 
^and/yoftheEGS. 



star name 


n/mas 


crjmas 


yu„/mas yr ' 


yujmas yr ' 


cr^„/mas yr ' 


cr^,,/masyr ' 


tr^/mas 


cr,, /mas 


RE J 0317-853 


34.380 


0.260 


-91.165 


-15.344 


0.435 


0.451 


0.3427 


0.2085 


LB 9802 


33.279 


0.238 


-78.894 


-27.041 


0.424 


0.412 


0.3042 


0.2030 


Refl 


4.62 


0.39 


10.76 


19.50 


0.782 


0.731 


0.4544 


0.1541 


Ref6 


3.51 


0.40 


0.00 


0.00 


0.000 


0.000 


0.4862 


0.4936 


Ref7 


1.57 


0.00 


-21.01 


-8.01 


0.698 


0.702 


0.4552 


0.3535 


Ref8 


0.17 


0.00 


0.00 


0.00 


0.000 


0.000 


0.4713 


0.3859 


Ref9 


1.23 


0.00 


0.00 


0.00 


0.000 


0.000 


0.4471 


0.2317 



ence field stars. At each of the three epochs, two HST orbits 
were used. Within each orbit, EGS Ir sequentially observed each 
star several times in a round-robin fashion for approximately 30 
seconds. The standard EGS data reduction algorithms ( |Nelan| 
& Makidonl|2002| l were employed to remove instrumental and 
spacecraft artifacts (such as photon shot noise, spacecraft jitter 
and drift, optical distortion of the EGS, differential velocity aber- 
ration, etc). The calibrated relative positions of the stars in each 
of the six visits were combined using a six parameter overlap- 
ping plate technique that solves for the parallax and proper mo- 
tion of ea ch star. This p rocess employed the least squares model 
GaussEit ( Jefferys et al.|1988 ) to find the minimum ;^^^ best solu- 
tion. 

The results of the EGS measurements for RE J 0317-853, 
LB 9802, and the reference stars are given in TablelS] The cr^ and 
cr,j are the 1 cr errors of the fit of the stars onto the ' master plate". 
Likewise, the parallax and proper motion errors are the 1 cr dis- 
persion in those values measured for the individual observations 



(e.g. RE J 0317-853 and LB 9802 were observed approximately 
four to five times each in every HST orbit, for a total of 24 to 
30 individual measurements). The errors quoted in TablelS] are 
typical of the performance of the EGS Ir (for comparison see 
e.g. Benedict et al.|[2007 l, indicating that our observations are 
nominal. The best solution was obtained by directly solving for 
the trigonometric parallax of Refl and Ref6, for which we ob- 
tain values consistent with their predicted spectroscopic paral- 
laxes. Likewise, we derive the optimal solution when we use the 
EGS Ir data to solve for the proper motion of Ref6 and Ref7. The 
bulk proper motion of the field is constrained by setting Ref6, 
Ref8, and Ref9 to have no proper motion. The astrometric ref- 
erence star Ref2 was not used because EGS Ir resolved it to be 
a wide binary system, which caused an acquisition failure in the 
second epoch. 

The parallaxes for RE J 0317-853 and LB 9802 differ by 
1.101 mas, which is about four times the 1 cr of their individ- 
ual errors. This result includes an application of the standard 
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fbf«)iKnjtanfrMiK K]^4i 




Fig. 2. The field of the binary WDs RE J 03 17-853 and LB 9802 
and the reference stars Refl, . . . , Ref9. Ref4 and Ref5 were later 
omitted since they were too faint. 



"lateral colour" correction that removes the apparent shift of 
an object's position in the FGS field of view due to the refrac- 
tive elements in the instrument's optical train. The correction is 
given as 6x = (B - V) ■ lex and 5y - {B - V) ■ ley. The co- 
efficients lex = -1.09 mas and ley - -0.68 mas are derived 
from the observed relative positions of two calibration stars, 
LATCOL^ (B - y = 1.9), and LATCOL_B (B - V = 0.2), 
at several HST roll angles. However, RE J 0317-853 is signifi- 
cantly hotter and bluer than the blue calibration star LATCOL_B. 
It is clear from Figure[T] that an error in the lateral colour cor- 
rection (especially in this case, along the FGS X-axis, which is 
nearly aligned with the line connecting the two circles marking 
the dates of the observations) will produce an error in the mea- 
sured parallax. To evaluate the validity of applying the standard 
lateral colour correction (which is based solely on a star's value 
of B-V) to RE J 0317-853, we revisited the interpretation of the 
astrometric results of the lateral colour calibration observations. 
Details of this "plausibility" investigation will be published as an 
STScI FGS Instrument Scientist Report (Nelan, in preparation) 
but summarized here. 

The spectral energy distribution (SED) of the two lat- 
eral colour calibrations stars, in addition to LB 9802, and 
REJ 0317-853 were convolved with the wavelength-dependent 
sensitivity of the FGS over its bandpass (the sensitivity decreases 
from ^ 20% at 4000A to ^ 2% at 7000A in a near linear 
fashion, where sensitivity refers to the probability that a pho- 
ton will be detected). The number of photons observed (i.e., 
actually detected) by the FGS for each star at a given wave- 
length {N photon si^)) was normalized to unity at (for the mo- 
ment) an arbitrary A,,. The effective FGS colour of each star 
was then defined to be the ratio of the wavelength weighted sum 
(TjU'^o - /I) * Nphotons{X)) for all A < Ao (the "blue" sum) to the 
similar "red" sum {YSiA- Ao)*N photonsi'^)) for all A > Ag over the 
FGS bandpass. The value of Ao is the boundary between the blue 
and red such that for a source emitting the same number of pho- 
tons at every wavelength the blue and red wavelength weighted 
sums are equal and the colour ratio is unity. For the FGS, we find 
that Ao = 5092A. 

The SEDs of both the red calibration star LATCOL^ and 
REJ 0317-853 were represented as black body curves with T - 
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Fig. 3. The spectrophotometric determinations of the absolute 
magnitude of the reference stars. The abscissa denotes the spec- 
tral type encoded in a way that FO corresponds to 0.0, GO to 
1.0, KO to 2.0 etc. and the spectral type subdivisions being given 
by the first decimal. The (red) open squares are the loci of main- 
sequence stars in this HR-diagram, and the (blue) open hexagons 
represent the giants (luminosity class III); the two curves show 
the resulting fits for both luminosity classes. The asterisks show 
the reference stars of this program. The spectral (sub)type was 
determined using low resolution spectra and the absolute mag- 
nitude was calculated using the fitted polynomial. 



2 900K and T = 50 000 K, respectively, while LATCOLJB and 
LB 9802 were represented by stellar model atmospheres using 
a code based upon the Kurutz models. For LATCOL_B, a so- 
lar abundance, Teff = 8 000K, and \ogg = 4.1 were assumed. 
For LB 9802, we assumed a hydrogen-atmosphere white dwarf 
with Teff = 16 030 K and log^ = 8.2. Using these SEDs, 
we computed for each star the wavelength weighted blue/red 
ratio described above, for which we found (blue/red) = 0.13, 
1.42, 1.79, and 2.54 for LATCOL_A, LATCOL_B, LB 9802, and 
REJ 0317-853, respectively. (A more precise estimate of the 
blue/red ratios for these four stars will use observed SEDs, which 
are currently unavailable. Here we simply evaluate the plausibil- 
ity of this concept.) 

If we assume that the lateral colour shift in the rela- 
tive position of two stars is proportional to the difference in 
their blue/red ratios, we can use the the astrometric results of 
the lateral colour calibration, which found that the blue star 
LATCOL_B was shifted by -1.87 mas relative to LATCOL_A, 
and their blue/red ratios to determine the proportionality con- 
stant a - -1.85/(1.42 - 0.13) - -1.44 mas. Applying this to 
REJ 0317-853 and LB 9802, we find the lateral colour-induced 
shift in the position of REJ 0317-853 relative to LB 9802 to 
be -1.08 mas. The parallax result cited in TablelS] already in- 
cludes a lateral colour correction of -0.25 mas in the position of 
REJ 0317-853 relative to LB 9802 (based solely upon the (B-V) 
of each star). This differs by -0.83 mas when using the difference 
in their blue/red ratios. If we apply this correction, the parallax 
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difference of the two stars is reduced to 0.27 mas, which is s; 1 cr 
of the individual measurements. We conclude that the two white 
dwarfs have the same parallax, and that this 0.27 mas difference 
is caused by the imprecise model SEDs used to construct the 
blue/red ratios. 

The measured parallax of LB 9802 is also affected by er- 
rors in the lateral colour correction, but to a lesser extent since 
at (B-V)= 0.07 it is closer to the colour of the blue calibra- 
tion star LATCOL_B (B-V= 0.2). Nonetheless, the parallax 
quoted in TablejS] may be too large by up to 0.4 mas, based on 
the difference in the predicted relative shift between two stars 
with (B-V) = 0.2 and 0.07 using the standard lateral correction 
and the blue/red ratio correction. Given the imprecision of the 
(blue/red) based correction, we take the parallax of LB 9802 to 
be ;7r = 33.279 ± 0.238 mas using the standard lateral colour 
correction. 

LB 9802 is 7" distant from RE J 0317-853 at a position an- 
gle P.A. = 145.856° as measured by FGS Ir. From the mea- 
sured proper motions (Table[5|, LB 9802 is moving away from 
RE J 0317-853 at 16.26 + 0.86 mas yr"' along a position angle 
of 133.62°, which is nearly aligned with the line of sight be- 
tween the two stars. The computation of the proper motions is 
dominated by the observations from the first and third epochs, 
which were performed at the same HST orientation. Therefore 
the uncertainty in the lateral colour correction has no effect. At a 
distance of 30.05 pc (calculated from the parallax of LB 9802), 
this corresponds to 0.489 ± 0.026AU yr^', i.e. 2.33 ± 0.12km 
s"'. We compare this tangential space velocity with an estimated 
orbital speed. If we assume that this is a bound binary system 
with a separation of 7" (210 AU at 30.05 pc), and with the total 
mass ranging from 2.02-2.31 Mq (see Sect. 4.1 1, a circular or- 
bit yields a period of 2004 yr (for the higher mass estimate) to 
2143 yr (for the lower mass); this corresponds to orbital speeds 
of 3.12-2.92 km s^^ for LB 9802 with respect to RE J 0317-853. 
These estimates are comparable to the tangential space veloc- 
ity measured by FGS. This result and the close spatial proxim- 
ity of the two stars supports the conclusion that LB 9802 and 
RE J 0317-853 constitute a bound system. 

Although the FGS photometry shows a peak-to-peak ampli- 
tude variation between V - 14.60 and V - 14.84 (with 0.01 ac- 
curacy estimated using LB 9802 as a reference) consistent with 
the result from Barstow et al. ( 1995| l, the sampling was not good 
enough to confirm the 725 second photometric variability quan- 
titatively by means of a Fourier analysis of this sparse data set. 

4. Determination of the stellar parameters 

4.1. Mass and radius determinations of RE J 031 7-853 and 
LB 9802 

To determine the mass of REJ 0317-853, we used synthetic 
bolometric colours and absolute magnitudes for carbon-oxygen 
(CO) core white-dwa rf cooling models with thick hydrogen la y- 
ers (Mh/M. = IQ-"*) ( |Wood|l995||Holberg & Bergeron|2006f 
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Fig. 4. Contour plots for \My - Mj^'=°|/mag as a function of 
mass in M0 and Tgff for CO (top), ONe (bottom) core com- 
positions constructed according to Eq.fT] and theoretical models 
from Wood ( 1995), Holberg & Bergeron (2006 ) for the CO mod- 
els, and Althaus et al. (2005 , 2007| for the"ONe models. The bai- 
to the right indicates the colour coding for the magnitude differ- 
ences, the line in the darkest region \M'^^ - My''"] < 0.5 mag 
delinating M"*"^ - My''° - and the vertical lines the possible 
range of effective temperatures (30 000-50 000 K). 



were calculated. The theoretical absolute visual magnitude was 
defined by 



Mf°(reff,m) 



-■ Mboi(T<,ff, m) - B.C.(reff , m). 



(1) 



when required, we used oxygen-neon (ONe) core white-dwarf 
cooli ng models with hydrogen layers of Mu/Mt = 10"^ (Althaus 
letal. 2005 2007 p] 

We determined the "observed" absolute visual magnitude 
M*^ = y-h51og;r-5 = 12.51 mag from V = 14.90 and 
TT = 0.033279". For a given effective temperature and surface 
gravity, the theoretical bolometric magnitude Mboi, the bolomet- 
ric correction B.C.=Mboi - My, and mass m for REJ 0317-853 

' http : //www . astro . nmontreal . ca/~bergeron/CoolingModels 
^ http : //www . f caglp . unlp . edu . ar/evolgroup/tracks . html 



The contour plots for \M^'^ - My''°\ are shown in Figured 
for the two possible core compositions. For both compositions, 
a satisfactory minimum could be reached only for parts of the 
range of effective temperatures between 30 000 and 50 000 K be- 
cause the tables were limited to an upper value of log g = 9.5 for 
the case of the CO cores (logg - 9.5 corresponds to a mass of 
1.37 M© for 30 000 K and a mass of 1.46 M© for 50 000 K) and 
to an upper limit of 1.28 M© for the ONe models. 

We calculated the minimum of IM*** - My''°\ for a given 
mass of our range of effective temperatures; when a mass solu- 
tion could not be reached inside the calculated grids, we extrap- 
olated the theoretical magnitudes. 
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Table 6. Mass and age estimations for RE J 0317-853 using dif- 
ferent core compositions and temperatures. The differences in 
radius estimates are caused by the different hydrogen content for 
different core models (see Sec. 14. 11). 



1.20 1.25 

M/'M 



Fig.5. The mass of REJ 0317-853 versus absolute V magni- 
tudes for an ONe white dwarf. The different curves correspond 
to the effective temperatures 30 000-50 000 K. Above 1.28 M©, 
we have to perform an extrapolation for T^-ff > 30 000 K. Since 
we cannot strictly estimate the extrapolation error, we visually 
added some uncertainty to the extrapolated values, which was 
subsequently used to estimate the errors in Tablel6] The red line 
denotes the "observed" My- 



For an effective temperature of 30 000 K, we estimated 
masses of 1.32 ± 0.02 (CO core) and 1.28 + 0.02 (ONe core). 
Our CO-core calculations are consistent with the estimates of 
1.31-1.37 M0 (Ferrario et al. 1997| l, who assumed a distance 
of 30 pc. The highest temperature for which we could obtain a 
solution in the IM"*"' - M*™| diagram is about 48 000 K from 
which we inferred a mass of 1.46 M©. Any additional extrapo- 
lation may introduce substantial uncertainty because we are then 
approaching the Chandrasekhar limit. 

In the grid of theoretical values for ONe cores, we performed 
significant extrapolation to obtain solutions above 30 000 K (see 
FigurelSll. For Teff - 30 000 K, we obtained a mass of 1.28 M© 
and interred an error of ±0.015 from the uncertainty in the ob- 
served visual magnitude and the parallax. For an effective tem- 
perature of 50 000 K, we derived 1.38 Mq with a slightly higher 
error estimate of 0.020 Mq due to the uncertainty of the extrap- 
olation. The results are summarised in Table|6] 

We applied the same procedure to LB 9802 by using our new 
parallax measurements and the information from the literature 
outlined in T ablefT] Our mas s estimate for the visual magni- 
tude given by Barstow et al. (|1995[l is consistent w ith the for- 
mer results ( [Ferrario et al.|1997||Kawka et al.|2007| see Table|7]l 
although we find that our calculations with the visual magni- 
tude provided by Kawka et al. ( 2007) is incompatible with our 
mass determination if we assume that the spectroscopicaUy de- 
termined masses for LB 9802 are correct. 

With the knowledge of the Mboi for a given mass, the radius 
can be directly estimated at a given Tgs. The radius estimates 
yield slightly different values when two core models are consid- 
ered (see Table. l6]|. This is caused assuming the assumption for 
the hydrogen layer mass to be Mh/M, = 10"'* in the CO cooling 
models ( Wood 1995) versus the Mh/M, - 10"^ content in the 
ONe cooling models ( Althaus et al.|2005 1. This produces differ- 
ent luminosities for a given effective temperature. 
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4.2. Age determination of REJ 031 7-853 and LB 9802 

The assessment of the cooling ages of REJ 0317-853 and 
LB 9802 is important to the understanding of the evolutionary 
history of the system. It was possible to evaluate the cooling 
ages of both objects with the mass estimates that we determined. 

For our estimations, we used the grids of white dwarf coohng 
sequences for CO and ONe cores (e.g. Wood|1995 Benvenuto 
& Althaus]! 1999 1 for their respective range of grid parameters; 
for masses above the available values, we extrapolated the age 
values in a way similar to that for the visual magnitudes (see 
Figs.|5]and|6]). 

Surprisingly, the difference in the cooling age of the two bi- 
nary components is smaller than formerly estimated. For both 
assumed chemical compositions, the cooling age of the non- 
magnetic white dwarf LB 9802 is within the error bars of the 
cooling age of the magnetic and very massive REJ 0317-853 
(see Tables[6]and[7]). For the case of an ONe core with an effec- 
tive temperature as high as 50 000 K, our conclusion is poorly 
constrained due to the extremely large uncertainties introduced 
by the extrapolation. 

Previous age estimates were unreliable because they inferred 
a coohng age of REJ 0317-853 shorter than that of LB 9802, 
simply based on its higher effective temperature. If we use the 
elementary theory of cooling by 'Mestel ( 1965 1 assuming for a 
fixed effective temperature of the white dwarf, the cooling age is 
a function of the mass and radius fcooi '^ M/R^. This means that 
the cooling age for low-mass white dwarfs (< 0.5 Mq) is simply 
proportional to mass M^^^. As the mass of the white dwarf ap- 
proaches the Chandrasekhar limit the radius asymptotically ap- 
proaches zero, which means that ages for a given effective tem- 
perature depend even more strongly on the mass. 

The masses estimated here are quite close to the 
Chandrasekhar limit (> 1.30 Mq) where post-Newtonian 
corrections should be considered for the stellar equilib- 
rium (Chandrasekhar 1964 Chandrasekhar & Tooper 1964[ l. 
However, these corrections mostly affect the dynamical sta- 
bility of the star, leading to collapse before reaching the 
Chandrasekhar limit, but induce only small corrections to mass- 
radius relationship. This is because the estimated radii are 
three orders of magnitude larger than the Schwarzschild-radius: 
GM/c^RwD ~ 10"^. Hence, we do not expect any effect on our 
mass determinations, as also noted by Koester & Chanmugam 
([1990}. 



5. The evolutionary history of the LB9802and 
REJ 0317-853 system 

The projected distance of 210 AU between the two white dwarfs 
and their small relative proper motion suggest that they are com- 
panions and therefore share a common origin. The ages of both 
objects should therefore be equal or comparable within the error 
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Table 7. Mass and age estimations for LB 9802 using different V 
magnitudes in the literature and an average effective temperature 
of 16000 K. 



V/mag mass/M0 fcooimg/Myr 



14.11' 
13.90^ 



0.84: 
0.76: 



0.05 
0.05 



279 
223 



lit 



' using the visual magnitude from Barstow et al. jl995f ; ^ using the 
visual magnitude frompKawka et al.,(,2007| 



bars; this condition must be fulfilled for the correct evolutionary 
schemes of both white dwarfs. 

The case of LB 9802 is straightforward because its evolution- 
ary history is not complicated by either a strong magnetic field 
or an extreme mass. Therefore, the simple single-star evolution 
of LB 9802 places constraints on the total age of RE J 0317-853. 

As mentioned above, previous analyses suggested a younger 
age for RE J 0317-853 than LB 9802 and for this reason the sys- 
tem was assumed to have an "age dilemma" (|Ferrario et al.| 
|1997| l. Therefore, an alternative scenario was proposed in which 
REJ 0317-853 is a result of the merging of two white dwarfs 
that have lower-mass progenitors. 

5.1. Single-Star origin of REJ 031 7-853 

With our new results, we undertook a more precise investigation. 
We firstly considered the single-star scenario for RE J 03 17-853. 
To determine the total age of LB 9802 and REJ 0317-853 from 
the zero-age main-sequence (ZAMS) to their current stage, 
we used the latest semi-empirical initial-to-final-mass relations 
(IFMR) HCasewell et"aL]pU09) [Salaris et al.1|2009l l to estimate 
their initial masses. By considering a diverse range of the the- 
oretical schemes to calculateg the IFMR (metallicity, overshoot 
parameter, etc.), we deduced the progenitor mass of LB 9802 to 
be in the range 4.0 - 4.5 M©. 

For the extremely high (final) mass of REJ 0317-853, the 
corresponding IFMR is quite uncertain. Theoretically, it was 
shown that 9-10 M© mass stars would evolve into massive 
oxygen-neon (ONe) white dwarfs because of the off-centred car- 
bon ignition in the partially degenerate conditions of their cores 
( |Ritossa et al.|1996HGa-cia-Berro et al.|1997| l. With these con- 
straints in mind, we consider more carefully a possible range of 
initial masses between 8 and 10 solar masses. 

The total age (time on the main-sequence plus the white 
dwarf cooling time) of LB 9802 depends strongly on its ini- 
tial mass. For the 0.84 Mq mass LB 9802, the initial masses in 
the range 4.0 - 4.5 M© yield main-sequence lifetimes of 170- 
130Myr (the progenitor ages were calculated using the evolu- 
tionary tracks from [Bertelli et alT (2009i for solar metallicity). 
This means that the total evolutionary age of LB 9802 is in the 
range 410- 450 Myr 

With 40-30 Myr, the pre-white-dwarf lifetime is extremely 
short for progenitor masses in the range between 8 and 10 Mq, 
respectively. For an effective temperature of 30 000 K and our re- 
sulting mass of 1.28-1 .32 M© for REJ 0317-853 (w hich would 
be the progeny of a 8 Mq star, see Casewell et al.|20 09; Salaris 
et al.|2009j ), we derive total ages in the range 320 - 340 Myr 

If alternatively we assume an effective temperature of 
50 000 K for REJ 0317-853 and a CO core, we end up with a 
total lifetime of a; 350 Myr; for the ONe core case, the corre- 
sponding value would be « 220 Myr We reiterate that our esti- 
mate for the errors is rather large in the ONe case at 50 000 K 



(see Table|6| because of the uncertainties in the extrapolation. 
Hence, omitting the case with ONe core at 50 000 K, we can say 
that the total age of RE J 0317-853 is in the range 320 -350 Myr. 
There are additional theoretical uncertainties in the IFMR 
due to magnetism and rapid rotation that should be important 
for an extreme case such as REJ 0317-853. The effect of both 
of these factors on the IFMR has been the subject of some dis- 
cussion. Dominguez et al. (1996) argued that rapid rotation has 



a positive effect on the core growth, such that a rapidly rotating 
star of mass 6.5 Mq may produce a white dwarf of mass 1.1- 



1.4 Mq. Observational evidence of this was found by Catalan 
'eraD('2008). REJ 0317-853 is the fastest rotating isolated white 
dwarf and this rotation may be a relic of a rapidly rotating pro- 
genitor. 

The assumption of a 6.5 Mq mass star as the progenitor does 
not relieve the "age dilemma" considerably since the progenitor 
age for this case is ~ 70 Myr, which does not differ much from 



the 40-30 Myr estimated for 8-10 Mq mass stars. Catalan et al 
( |2008| l also argued that MWDs are relatively more massive than 
expected on the basis of their inferred progenitors via the IFMR 
of non-magnetic white dwarfs. However, |Wickramasinghe & 



Ferrario ( 2005 ) and [Ferrario et al.| ( |2005| l both concluded based 
on their population synthesis studies this effect is of only minor 
importance. Since the effect of rotation and magnetism on the 
evolutionary age is unclear or rather small, we did not consider 
them in our age estimations. 

Based on these considerations, we conclude that the total age 
of LB 9802 is 410 - 450 Myr at least ~ 100 Myr older than the 
respective value for REJ 0317-853 (320 - 350 Myr). This dis- 
crepancy implies that the single-star evolution scenario might 
not be applicable to REJ 0317-853. 

However, the mass estimates leading to the cooling ages de- 
termined above neglected the influence of magnetism. The mag- 
netic nature of REJ 0317-853 is likely to affect the determina- 
tion of its mass because of the mass-radius relation. lOstriker &1 
Hartwick^(1968) discussed the effect of magnetism and rapid ro- 
tation on white dwarfs. Both magnetism and rotation act against 
the gravitation, causing an extended radius; hence, white dwarfs 
with strong internal magnetization have larger radii for a given 
mass. 

To calculate the cooling tracks from synthetic colours and 
magnitudes of white dwarfs, mass-radius determinations are 
used implicitly. Hence our estimates of the masses and ages are 
impaired by the lack of mass-radius relations taking into ac- 
count the effect of the magnetic field. For a white dwarf with 
1.05 Mq, the radius is increased by a factor e^^ = e-''^*, where 
6 is the ratio of the magnetic energy to the gravitational energy 



of the star, and n is the polytropic index ( Shapiro & Teukolsky 



,1983) . In the case of REJ 0317-853, an internal magnetization 
of < B >= 10^^ - 10^^ G seems plausible; this would imply that 
6 ~ 0.1 and therefore an increase in the radius by ~ 40%. Since 
REJ 0317-853 has an even higher mass, n is in this case close 
to 3 and thereby the increase in radius for a given mass is even 
higher 

For an effective temperature of 30 000 K, our measured ra- 
dius is 0.410 X 10^2 Rq^ whereas for 50 000 K it is 0.295 x 
10"^ Rq. When we correct the influence of the magnetic field 
on the radius we end up with a higher mass than determined in 
Sect. 4.1 If RE J 0317-853 were of higher mass the cooling time 



would increase so that the age dilemma no longer exists for the 
assumption of single-star evolution. 

As an initial consideration, cooling ages of ~ 400 Myr, 
which would diminish the age inconsistency, are possible for 
REJ0317-853, ifithasamassof 1.32 Mq rather than 1 .28 Mq 
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(ONe case; 0.04 M0 discrepancy), or 1.38 M0 rather than 
1.32 Mq (CO case; 0.06 Mq discrepancy), for an effective tem- 
perature of 30 000 K. The corrected radius of Rq - 0.32 x 
10"^ R© implies a mass of 1.38 M© from the mass-radius re- 
lationship. This value implies that the corrections are plausibly 
high enough to account for the missing evolutionary age as dis- 
cussed above. 

If we consider T^s = 50000K for REJ 0317-853, the mass 
estimates based purely on the total evolutionary age of the sys- 
tem would imply values well above the Chandrasekhar hmit. 
Although it is known that strong internal magnetic field strengths 
also modify the Chandrasekhar limit (Ostriker & Hartwick 



|1968) , it is still difficult to quantitatively assess the masses and 
their eff'ect on cooling ages in this regime. 



5.2. Binary origin of REJ 031 7-853 

The merger scenario for ultramassive white dwarfs was initially 
propos ed by Bergeron et al. ( 1991 1 for GD50, Marsh et al. 
( |199 7) proposed that this scenario could explain the properties of 
the hot and massive white dwarf population. For RE J 0317-853, 
it was similarly proposed to explain both the high angular mo- 
mentum and high mass of this star ( Ferrario et al.|1997 1. Vennes 
et al.| ([2003 ) also suggested that the scenario could produce 
a strong and non-dipolar magnetic field. They argued qualita- 
tively that the high angular momentum is a result of the total or- 
bital momentum of a coalescing binary and that the strong non- 
dipolar magnetic field can be generated by dynamo processes 
due to the differential rotation caused in turn by the merging. 

The type of binary evolution that can lead to a double- 
degenerate system has been investigated in detail, since it repre- 
sents a channel for producing SN la explosions (|Webbink|1984 



Iben & Tutukov|1984i. In this scenario, a binary system consist 



ing of two intermediate-mass stars (5-9 Mq) goes through one 
or two phases of a common envelope (CE) and evolves to a dou- 
ble white dwarf system. If the final double-degenerate system 
has orbital periods in the range between 10 s and 10 h, it will lose 
angular momentum through gravitational radiation and merge in 
less than a Hubble time. The merging process leads to a massive 
central product with a surrounding Keplerian disk. Depending 
on the total mass of the system, the temperature in the envelope 
and the accretion to the merger product, the system can evolve 
either to a SN la or by an accretion-induced collapse (AIC) to a 
neutron star. When the total mass of the system is insufficient to 
create the density and the temperature to burn carbon under de- 
generate conditions, the system will end up as an ultra-massive 
white dwarf. 

To test whether this scenario is indeed applicable to the case 
of RE J 0317-853, we have to trace back to the point in the stel- 
lar evolution where the merging could have happened, using the 
cooling age of REJ 0317-853 and subtracting it from the total 
evolutionary age of LB 9802. Using this progenitor age estimate 
and the theoretical constraints from the theory of binary star evo- 
lution, we can estimate the masses of the possible merging coun- 
terparts. 

We begin by estimating the mass of the (secondary) binary 
component that needs longer to become a white dwarf. To ob- 
tain a lower limit to its mass, we assume the longest time from 
the main-sequence to the merging process considering the mass 
transfer episodes predicted by the binary scenario. After both 
white dwarfs are formed, the time needed for the binary to merge 
due to gravitational radiation depends strongly on the orbital pa- 
rameters and mass of the double-degenerate system. Depending 
on the properties of the system, coalescence can be as fast as 



0.1 Myr or as slow as 200 Myr ( |Iben & Tutukov|1984| l. To ob- 
tain a lower limit to the total evolutionary time for the system, 
we neglect the time needed for the double-degenerate system to 
coalesce. 

|Iben & Tutukov] ( |1985[ ) discussed the evolution of 3 to 
12 Mq stars that experience two phases of mass transfer The 
phase of the mass transfer can take as long or even longer than 
the time the star spends on the main-sequence. For a 5 Mq star, 
the main-sequence phase lasts ~90Myr ( Bertelli et al.|[2009) , 
while in the binary-evolution scenario it takes 140 Myr from 
the main sequence until the formation of the white dwarf. This 
means that 230 Myr are needed for a 5 Mq star to evolve into a 
white dwarf rather than the 100 Myr that we assumed for single- 
star evolution. 

The possible cooling ages considered for LB 9802 (280 Myr) 
and RE J 0317-853 (280- 320 Myr, when we assume aneff'ective 
temperature of about 30 000 K) imply that the maximum time 
needed for binary evolution is at most the main-sequence age 
of LB 9802, which is 130-170 Myr (for 4.0-4.5 Mq). The upper 
limit of 170 Myr is comparably short relative to the 230 Myr of 
binary evolution time. This provides a lower mass limit for the 
system. The resulting mass of a white dwarf that is a product of 
a 5 Mq star in this binary evolution scheme is 0.752 Mq (Tbenl 
|& Tutukov|1985[ ), which is lighter than inferred from the IFMRs 
determined for single-star evolution. 

Since the pre-white-dwarf evolution is too long for an initial 
5 Mq, star we need a more massive progenitor hence should end 
up with a secondary white dwarf more massive than 0.752 Mq. 
For the primary star, we assume that it has only a slightly higher 
mass than the secondary to deduce a lower limit to the total coa- 
lescing mass. 

However, this assumption leads to serious inconsisten- 
cies, because the total mass of two components would result 
in more than 1.5 Mq being above the masses estimated for 
REJ 0317-853. This lower limit is also robust when we consider 
mass loss. Firstly, smoothed particle hydrodynamic (SPH) sim- 
ulations show that only a very sma ll mass loss is expected dur - 
ing merging (~ 10"^ Mq, see e.g. [Loren-Aguilar et al.|[2009 1, 
and secondly, we expect almost all of the Keplerian disk to 
be accreted on the merger product. Wind mass-loss from the 
Keplerian disk is assumed to be lower than 10% of the accre- 
tion rate ( |Mochkovitch & Livio|| 1990) 1; this means that at least 
90% of the disk is expected to be accreted. "Loren-Aguilar et al. 



( 2009| also estimate 0.1 - 0.3 Mq for the disk masses, which 
would imply a total mass loss of < 0.01 - 0.03 Mq. 

We note that infrared studies have detected possible disks 
surrounding massive white dwarfs ( Hansen et al.| |2006). This 
included REJ 0317-853 for which no convincing evidence of 



a disk was found in the Spitzer IRAC bands (see also |Farihi 
2008) l. If REJ 0317-853 were the product of a merger 



et al. 



of two white dwarfs, all of the matter from the Keplerian disk 
should have been accreted. In this scenario, total mass limits 
well above the estimated RE J 03 17-853 mass cannot be avoided. 
This estimation eliminates the possibility of a binary origin for 
REJ 0317-853 with a current effective temperature as low as 
30 000 K. 

However, if the total mass of the binary system does not ex- 
ceed the estimated value for RE J 0317-853, the time needed for 
the accretion of all the material from the disk is much longer than 
the evolutionary timescale. The accretion rate is expected to be 
< 10"'^ MQ/yr for flows with laminar viscosity (Loren-Aguilar 



et al.|2009 ). For disk material of 0.1 - 0.3 Mq, that its complete 
accretion time of 1 - 3 x 10^ Myr is three orders of magnitude 
longer than the evolutionary timescale. 
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If the binary scenario were correct, the Keplerian disk should 
have been observed unless the accretion rate of the disk was 
much higher than theoretically predicted. Only accretion rates 
higher than 10"'" M0/yr would lead to a total disappearance of 
the disk. 

When two equal-mass white dwarfs merge, the symme- 
try of the process leads to a rotating ellipsoidal composed of 
CO around the white dwarf rather than a Keplerian disk. If 
RE J 0317-853 were still in the process of accretion we would 
have observed CO in the spectra but this is not the case. On 
the other hand, if all the material of the surrounding ellipsoid 
had already been accreted (mass-loss can be neglected as dis- 
cussed above) the mass of RE J 0317-853 would be higher than 
observed (above the Chandrasekhar limit). 

We also considered the possible effect on the cooling ages 
of additional heating of the white dwarf core due to the merging 
process. Recent SPH simul ations indicate the possibihty of heat- 



ing to ~ 1() ^ K in the core ( Yoon & Langer 2005 Loren-Aguilar 

^'-^ dependence of the 



et al.[2009 1. However, because of the T~ 
cooling age according to the elementary theory (Mestel"1965i), 
the effect of this extra heating on the cooling ages is expected to 
be small (~ 2Myr) and can be neglected. 

When we consider the ONe core case for an effective tem- 
perature of 50 000 K, leading to an average cooling age of ~ 
190Myr, we end up with an upper limit to the evolution time 
of the secondary of 220-260 Myr; the 40 Myr spread in evolu- 
tionary time is only due to the uncertainty in the LB 9802 pro- 
genitor mass (between 4.0-4.5 Mq). Our estimated upper limit 
to the total age is comparable to the evolutionary timescale of a 
5 M0 star in a binary system as considered above. However, the 
cooling time estimate for REJ 0317-853 in this case is consid- 
erably uncertain (see Tablel6|l due to the extrapolation. Within 
these large error margins, we would in principle be able to ob- 
tain a sub-Chandrasekhar mass for the merger product, but this 
process is very unlikely when we consider the time needed for 
the white dwarfs to merge (10-100 Myr |Iben & Tutukov|1984| l. 
Nevertheless, the possibility of a binary origin for an ONe core 
RE J 0317-853 at Tetj = 50 000 K cannot be entirely excluded. 

We note that the effect of magnetic field strength on the struc- 



ture of the white dwarf, considered in Sect. 5. 1 is also important 



to binary evolution. The implementation of this effect leads to an 
inference of slower cooling for REJ 0317-853 as in the single- 
star scenario. This would yield shorter progenitor timescales for 
a constant evolutionary time, leading to the lower limits on the 
total mass of the coalescing white dwarfs becoming even more 
massive. This diminishes again the probability of binary evolu- 
tion for Teff = 30 000 K. However, for Teff = 50 000 K the uncer- 
tainties still permit the possibility of merging. Furthermore, the 
effect of magnetism on the stellar structure ensure that this sce- 
nario remains favourable due to the higher Chandrasekhar mass 



limit ( |Ostriker & Hartwick|1968| [Shapiro & Teukolsky|1983 1 



6. Discussion and conclusions 

REJ 0317-853 belongs to the very rare population of ultra- 
massive white dwarfs with masses exceeding 1.1 M0. The com- 
peting theoretical explanations of the origin of these white 
dwarfs are single-star evolution versus the merging of two de- 
generate stars. Without considering mass-loss during stellar evo- 
lution, we have shown that an upper limit of 1 . 1 M0 for the final 
white dwarf mass would exist for the white dwarfs because of the 
ignition of carbon in the core of the progenitor star. However, 
taking into account the effect of mass loss, high-mass ONe-core 




Fig. 6. The mass of REJ 0317-853 versus logarithmic age in 
years for an ONe core white dwarf. The different curves cor- 
respond to the effective temperatures 30 000-50 000 K. Since we 
cannot strictly estimate the extrapolation error we visually added 
some uncertainty to the extrapolated values, which was subse- 
quently used to estimate the errors in Tablel6] 



white dwarfs can be produced (see |Weidemann||2000 for a re- 
view). Furthermore, it was proposed that even 9 to 10 Mq mass 
stars evolve into ONe core white dwarfs of mass 1.26 and 1.15 
respectively, because of the off-centred carbon ignition in the 
partial degenerate conditions of their cores ( Ritossa et al.|1996 
I Garcia-Berro et al.|1997| i. 

In the light of our current results, we have undertaken a more 
precise investigation of the possible evolutionary scenarios for 
REJ 0317-853. We have shown that the cooling ages are almost 
the same for the two components. The detailed analysis very 
much depends on a precise determination of the effective tem- 
perature; for Teff - 30 000 K, we can use the calculations by 
Wood] ( |1995| l and [Benvenuto & Althaus| ( |1999 ) and conclude 
that within the limits of the uncertainties REJ 0317-853 is at 
least as old as LB 9802. For a consistent interpretation of the 
system, we also have to take into account the time scales of 
the pre-white-dwarf evolution. The more massive progenitor of 
REJ 0317-853 should evolve more rapidly than the progenitor 
of LB 9802. Taking this into account, the total age difference 
between LB 9802 and REJ 0317-853 amounts to ~ 100 Myr if 
single-star evolution is considered. 

On the other hand, the alternative binary merger scenario 
proposed by Ferrario et al. (1997) and Vennes et al. ( 2003 1 as 
a solution to this age dilemma has severe drawbacks. When 
the evolutionary timescales are considered, the progenitor age 
of REJ 0317-853 at T^s = 30000K yields lower Hmits on the 
mass of the merger product that is considerably higher than its 
estimated mass for all cases. For REJ 0317-853, we have large 
uncertainties in the cooling age estimate only for an effective 
temperature of 50 000 K, so that we cannot fully exclude the bi- 
nary scenario. 

We have also considered the effects of the magnetic fields on 
both of the scenarios. Magnetic fields cause an increase in radius, 
hence an underestimate of the mass, which would imply longer 
cooling ages than estimated. For the case of Teff = 30 000 K, 
the effect of magnetism makes the single-star scenario possible 
while further eliminating the binary merger origin; for the high 
Teif of 50 000 K, even the inclusion of magnetic effects ensures 
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that the single-star scenario is possible; the binary scenario re- 
mains possible within our large uncertainties. 

With our measurement of the parallaxes and relative proper 
motion of RE J 0317-853 and LB 9802 with HSTs FGS, we 
have established that the wide binary system of these two stars 
is indeed a bound system. We have estimated the masses and 
ages of RE J 0317-853 and LB 9802 based on the current white 
dwarf cooling tracks for different core compositions and hydro- 
gen layer masses. Owing to the magnetic nature of this object, 
the temperature determination of RE J 0317-853 is difficult and 
should be repeated in the future taking into account all available 
observations and including a more detailed determination of the 
magnetic field geometry. 

For the mass and radius determination, we have consid- 
ered the highest and lowest possible temperature and with these 
estimates we have discussed the evolutionary history and the 
possible origin of REJ 0317-853. Our results show that for a 
cooler, less massive REJ 0317-853 the binary scenario can be 
excluded within our uncertainties. We also proposed that the 
"age dilemma" might be solved when the effects of the mag- 
netism on the structure of the white dwarf is considered. If 
REJ 0317-853 were hotter and more massive, then a binary ori- 
gin scenario would be more plausible. 
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